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Abstract 

The aim of this paper is to present basic notions of dense mat- 
ter physics and some of its applications to geophysics and astron- 
omy. Topics covered in the paper include:basic observational data, fun- 
damental ideas of static high pressure experiments, notions of theoret- 
ical dense matter physics, and finally some details about theoretical 
work on dense matter physics and its astronomical applications in 
Serbia. 

Introduction 

Astronomy is a science full of contrasts. Nearly all physical quantities 
measurable in astronomy reach values which are hardly immaginable on 
Earth. Two interesting extremes which occur in space science concern the 
mass densities of the objects studied: they span about 30 orders of magni- 
tude. At the lower limit, astronomers encounter various kinds of nebulae, 
the interstellar and intergalactic medium, where the density can be as low 
as 10~ 15 kg m~ 3 . Near the upper limit are various kinds of objects in which 
matter is subdued to high values of pressure (and temperature) and where 
the density can go upwards to values of the order of 10 15 kg m~ 3 . 

The aim of this lecture is to review the basic notions of dense matter 
physics and its applications in astronomy. The field is huge, and the present 
lecture is far from being a complete review of the field. The topics which 
are covered are listed in the abstract; they were chosen so as to reflect the 
author's research experience but combined in a way to make this lecture an 
introduction to the basic notions and some recent research results. 
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Basic data 



Celestial objects in general, and those in the planetary system in par- 
ticular, are observable since prehistoric times. The nature of space science 
has dramatically changed since the invention of the telescope, as attempts 
to measure and draw conclusions steadily gained over pure contemplation of 
phenomena. 

The telescope, in combination with the laws of celestial mechanics, gave 
the possibility of determining masses and radii of various objects. Using these 
data, it became possible to calculate the mean densities of objects observed. 

Starting from the mass M and radius R of an object, and assuming spher- 
ical shape, it is clear that the volume is V = (4/3)7rr 3 , while its density is 
p = M/V . This kind of calculation gave the first possibility of gaining some 
knowledge on the chemical composition of other planets. Densities of various 
planets, derived from the observed data, could be associated with various ma- 
terials known on Earth. Modern values of planetary densities, taken from the 
National Solar System Data Center at http://nssdc.gsfc.nasa.gov/planetary 
are presented graphically in Fig.l. Planetary densities are mutually widely 
different: from those which are rarer than water to those which are 5-6 times 
denser. This is a consequence of a set of processes in which the planetary 
system was formed about 4.6 billion years ago. It is usually taken that the 
present distribution of planetary densities is a consequence of the spatial 
distribution of temperature and chemical composition of the material in the 
protoplanetary disk. For an example of recent work on planetary formation 
see P and references given there. 

Real experimental work in planetary science became possible in the sec- 
ond half of the XIX centurywith the developement of seismology. The 
physical idea on which seismology is based is extremely simple. A point in 
the interior of the Earth in which an earthquake occurs becomes a source 
of waves which propagate along the surface and through the interior of the 
planet. Measurements of the period, amplitude and damping rate of the os- 
cillations give information about the earthquake which is the source of the 
elastic waves, but also on the physical conditions along their path. 

Figure 2 is a historic example of a seismogram j2j. It shows the NS 
seismic waves recorded in Gottingen (Germany) during the earthquake in 
San Francisco on April 18,1906. The fact that the quake was recorded by the 
relatively crude equipement of the time shows the sensitivity of the method, 
but also the power released in this event. The symbols P and S denote the 



2 



6000 



>• 3000 - 



20 30 

mean distance (A.U.) 



Figure 1: Densities in the planetary system 

compressional and shear waves respectively, while PP and SS are waves which 
bounced 1 or 2 times from the surface. 

The main data which can be deduced from the analysis of seismograms 
are the values of the speed of propagation of different waves. Analyzing 
the wave profiles gives the possibility of determining the position of various 
discontinuities deep in the core [H|. The measured values of the speeds of 
propagation open the possibility of determining the compressibility of the 
mixture of materials through which the wave propagates. It is known that 
in a liquid the speeed of a compressional wave is related to the adiabatic 
compressibility K$ in a simple way: 

2 K S /.x 

v P = — (1) 

where 

The logic behind the applicability of these two expressions in planetology is 
simple: in the case of the Earth, it is usually assumed that its outer core is 
liquid. Accordinglyby measuring the speed of seismic waves resulting from an 
earthquake, eqs. (1) and (2) give the possibility of drawing some conclusions 
concerning the equation of state of the material. These expessions are based 
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Figure 2: NS waves in the San Francisco earthquake of 1906. 

on a simplified model; note that real work in seismology is based on more 
complicated equations, but the reasoning behind remains the same. 

Laboratory experiments: a few notes 

Historians of science have discovered that a certain wealthy English gen- 
tleman named Mr. Canton sometime in the XVIII century made a first at- 
tempt to perform an experiment under high static pressure. He compressed 
water at room temperature to a pressure of the order of 0.1 GPa. To his 
astonishement water was transformed into ice [IJ. We know today that the 
phase diagram of ice under high pressure has at least 12 phases jH]. 

Systematic experimental work on materials under high pressure has been 
initiated towards the end of the XIX century by P.W.Bridgman at Harvard 
University jH] . For his life- long work Bridgman was awarded the Nobel prize 
in physics in 1946. In his experiments, Bridgman used presses which could 
contain large samples and had small P — T gradients, which was very useful 
for the measurement process. At the same time, the accessible region of the 
P — T plane was extremely small, and the presses were expensive to build 
and complicated to operate. 
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A big breakthrough occured in the middle of the last century with the 
invention of the diamond anvil cell (DAC). Details about the construction of 
this instrument, the choice of diamonds and some of the experimental meth- 
ods are available in the literature such as [7j. An interesting step forward 
has occured a short time ago: the so called "ruby scale" which is used in 
nearly all measurements in DACs has been corrected, and the validity of this 
modification tested to P < 100 GPa,T < 850 K |H|. 

Examples of experimental studies 

Drawing conclusions about the behaviour under high pressure of astronom- 
ically interesting materials is comlicated by the experimental methodology 
and the complexity of phenomena occuring with the increase of the den- 
sity. An excellent example, spanning more than 70 years of research and still 
unsolved, is offered by the behaviour of hydrogen under high pressure. 

The problem of dense hydrogen can be resumed in the following way. 
It has been predicted numerous times in the last 70 years that hydrogen 
should become metallic at a pressure of the order of 250 — 300 GPa. These 
predictions were made by different authors and theoretical methods. As 
they all were giving the same order of magnitude, it was taken almost as an 
experimental result. However, when such values of pressure became acces- 
sible, and experiments performed [9j, hydrogen did not become a metal for 
P < 342 GPa. More recent experiments in a DAC have shown that molec- 
ular hydrogen continues to exist up to P = 316 GPa and that it becomes 
black at P = 320GPa ^U]. It was also experimentally obtained that the 
semiconductor —> metal transition in fluid hydrogen occurs at P = 140 
GPa; T = 3000 K |TT|. Finallynote that recent theoretical work indicates 
the possibility that hydrogen turns into a metal only at P > 400GPa [T2"] . 

In the field of shock compression experiments, the situation also seems 
to be stationary. At the time of this writing (October ,2005) there are no 
new results from the group in the Lawrence Livermore National Laboratory 
(LLNL) performing laser shock experiments. Some experimental work is 
going on in Rochester, but the next opportunity to obtain " big" results will 
have to wait for the completion of the National Ignition Facility (NIF). This 
facility is being built in LLNL, and it is expected to be completed in 2008. In 
the meantime " flying plate" experiments are being performed in Sandia 
Labs. (USA) and in Arzamas (Russia). They have reached 1.4 MBar,but 
their results for the equation of state (EOS) of hydrogen do not agree with 
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previous shock compression results from Livermore ^3] • Space limitations do 
not allow going into further details, but the point is that although hydrogen 
is the simplest chemical element, the most abundant in the Universe, its high 
pressure behaviour is extremely hard to predict. 

Another closely related material is water. Its importance for the existence 
of life on Earth is clear. In planetologyits importance stems from the fact 
that water ice (or ices containing water and different other materials) exist 
on some planets and their satellites and as such contributes to their albedos. 
This implies that knowledge of the phase diagram of ice(s) is necessary for the 
correct interpretation of the observed values of albedos of various planetary 
system objects. The "complexity" arises from the fact that the phase diagram 
of ice has at least 11 parts. The physics of ice is of considerable importance in 
studies of interstellar grains. See, for example. |14j and references given there). 

Studies of materials under high pressure are important for geophysics. 
Just as an example, note a recent experiment in which the ellasticity of hep 
iron has been studied in a DAC at room temperature for P = 112 GPa. The 
anisotropy of the velocity of sound has been measured, and the result is in 
close correspondence with the anisotropy of the Earth's inner core j!5j . 

A note on theoretical work 

Technological developement has made possible the enlargement of the 
experimentally accessible part of the P — T plane. However, an incomparably 
greater region of the same plane can be studied theoretically. Any specimen 
of a real material represents a typical example of a many body problem. In 
rigorous terminology of statistical physics, the many body problem is defined 
as the eigenproblem of the following Hamiltonian 

N *2 N N 

h = £(-£r) v? + E v{Si) + E - (3) 
t=i /m i=i i,j=i 

In this equation N denotes the number of particles in the system, m is 
their mass and all the other symbols have their usual meanings. In all real 
examples of many-body systems,the sums in this expression are impossible 
to calculate exactlybecause the number of particles N is of the order of 
10 23 . The art of theoretical dense matter physics in fact consists in finding 
methods of calculation of sums in eq. (3), because this ultimately leads to the 
possibility of determination of the thermodynamical potentials (and phase 
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transition points) of a system. Attempting to review the whole field would 
be out of place, so it the following we shall concetrate only on a contribution 
to dense matter physics by researchers in this country. 

The contribution from this country 

At the beginning of the sixties, Pavle Savic and Radivoje Kasanin have 
started developing a semiclassical theory of dense matter. Their theory is 
called " semiclassical" because it uses some notions of classical in combination 
with some ideas of atomic physics. The starting point of their research was 
trying to find a simple relationship between the solar and planetary mean 
densities, star ting form the observed data known at the time. It turned out 
that such a relationship exists and has the following form: 

P = P^ (4) 

where po is the mean solar density and <p is an integer. Choosing tp 6 [—2, 2] 
it becomes possible to reproduce the observed planetary densities. Although 
this fit "works" Savic and Kasanin did not propose a physical mechanism 
which could be at its origin. Using this result and a known fact from both 
geophysics and laboratory high pressure experiments that at some values of 
pressure abrupt changes of the density of materials occur,they qualitatively 
concluded that the electronic structure of materials changes under high pres- 
sure [To] . 

For the sake of historical completeness, it should be noted that they were 
not the first to propose such an idea;the first to advance the idea of "atomic 
destruction" under high external pressure was P.W.Bridgman back in 1927 
17j. A couple of years after, Fermi played with the same idea,while exploring 
various possibilities offered by the (then new) Schrodinger equation. 

Starting from these ideas, Savic and Kasanin proposed a set of 6 exper- 
imentally founded postulates which govern the behaviour of materials un- 
der high pressure. Developing these postulates, they have set up a calcula- 
tional scheme, which gives the possibility of theoretical studies of dense matter 
physics, both in laboratory and astronomical applications. Details are given 
in (for example) [18 . 

Input data needed for the modelling of the internal structure of a planet, 

satellite or asteroid within this theory are only the mass and the radius of 
the object. Starting from this couple of values, it becomes possible to calulate 
the values of the following parameters: 
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depth [km] 


0-39 


39-2900 


2900-4980 


4980-6371 


p max [kg m- 3 ] 


3000 


6000 


12000 


19740 


Pmax [GPa] 


25 


129 


289 


370 


Pmax [K] 


1300 


2700 


4100 


7000 



Table 1: The interior of the Earth 



• the number of layers in the interior and their thickness; 

• the distribution of P, p, T within the layers; 

• the mean molecular mass of the chemical mixture that the object is 
made of; 

• the magnetic moment of the object; 

• the limits of the physically allowed interval of values of the angular 
frequency of rotation of the object. 

All planets except Saturn and Pluto have been modelled. Models of the 
Moon, Galileian satellites, the satellites of Uranus, Triton and Titan, as well 
as the asteroids 1 Ceres and 10 Hygiea have been calculated so far. The 
first celestial body to be modelled within this theory was the Earth; as an 
illustration,table 1 contains the main results on the interior of our planet. 

The algorithm used in the theory proposed by Savic and Kasanin is ex- 
tremely simple. In spite of that, results of its applications to planetological 
problems are in good agreement with more complex models. A model of the 
crust, called CRUST5.1, used by the U.S. Geo logical Survey, is acessible at 
the following internet address http:/ /quake. wr.usgs.gov/research/structure/ 
CrustalStructure/index.html. In most parts of the Earth, values of the thick- 
ness of the crust are in good agreement with those given by the model of 
Savic and Kasanin. Disagreement is obvious near the poles,but finding its 
explanation is at present an open question. 

For all the planetary objects modelled within this theorythe value of 
the mean atomic mass A of the material they are made of was also calcu- 
lated. These values, collected from various papers, are presented in table 2. 

Interesting conclusions can be drawn from table 2. Triton and Neptune 
are not made up of the same material, which agrees with known results of 
celestial mechanics that Triton is a captured body. Differences of A are also 
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object 


mean atomic mass A 


satellite 


mean atomic mass A 


Sun 


1.4 


Moon 


71 


Mercury 


113 


Jl 


70 


Venus 


28.12 


TO 

J2 


71 


Earth 


26.56 


J3 


18 


Mars 


69 


J4 


19 


1 Ceres 


96 


Ul 


38 


Jupiter 


1.55 


U2 


43 


Saturn 


/ 


U3 


44 


Uranus 


6.5 


U4 


32 


Neptune 


7.26 


U5 


32 


Pluto 


/ 


Triton 


67 



Table 2: The composition of the planetary system 



visible for the Earth-Moon system,but their explanation is more complex. 
Modern work on the origin of the Moon indicates that it is a result of the 
colision of a body of the size of Mars with the Earth, which happened about 
4.5 billion of years ago. The necessary condition for the formation of the 
Moon is that the impact produces solid and liquid debris. Debris from this 
impact formed a disc around the Earth that coalesced to form the Moon. 
Recent numerical simulations jTH] show that the colision was slow,with the 
impactor moving with a speed less than about 15 kms^ 1 . According to 
these results, the deep interior of the Moon contains "the remnants" of the 
impactor, while the surface layers are in fact debris from the impact. Judging 
by their values of A, Mars and Triton appear to be chemically similar. Using 
only this similarity, one could immagine that Triton formed near the present 
orbit of Mars,was somehow ejected from this region and captured by Neptune. 
This idea demands "fine tuning" of the relative velocities of all the objects 
in the problem,which is difficult to achieve. The relative differences of the 
average atomic masses are visibly large, and they are certainly mostly caused 
by real physical differences between various objects but to some extent also 
by the extreme simplicity of the theory. 

A few words are in order on laboratory applications of this theory. A 
detailed algorithm has been developed for the calculation of the pressure at 
which a first order phase transition can be expected in a material under high 
pressure. It has been applied to about 20 materials, chosen at random,for 



9 



which reliable experimental data could easily be found in the literature |2*Uj . 
It was shown that discrepancies with experimental results exist, but that they 
vary between practically zero and 30 percent, and that they can be explained 
by differences in the structure of the materials studied and the nature of the 
inter and intramolecular potentials in them. These results are interesting on 
their own, and they also lend weight to all astronomical applications of this 
theory 

Conclusions 

Dense matter physics is a fascinating field of pure physics with interest- 
ing applications in space science.lt is practiced in this country for several 
decades, and the results are very encouraging. The aim of this contribution 
was to review the basic notions and some results of dense matter physics in 
general, and obtained in this country in particular. 
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